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Abstract 

We consider two closely related Riccati equations of constant parameters 
whose particular solutions are used to construct the corresponding class of 
supersymmetrically-coupled second-order differential equations. We solve an- 
alytically these parametric periodic problems along the positive real axis. Next, 
the analytically solved model is used as a case study for a powerful numerical 
approach that is employed here for the first time in the investigation of the en- 
ergy band structure of periodic not necessarily regular potentials. The approach 
is based on the well-known self-matching procedure of James (1949) and imple- 
ments the spectral parameter power series solutions introduced by Kravchenko 
(2008). We obtain additionally an efficient series representation of the Hill dis- 
criminant based on Kravchenko's series. 
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1 Introduction 

Natural periodic and quasiperiodic structures have drawn the attention of mankind 
since old times pQ. Nowadays, in the technological world, the scientists, espe- 
cially those involved in the area of material science, bring forth at a tremendous 
pace new artificial specimens in which finite periodic structures are the main 
components. A rich mathematical background related to periodicity has been 
developed along the years [21 [3]. The construction of new periodic potentials and 
the analysis of their specific properties could be valuable in guiding the mod- 
ern technological design. Within this context the relationship between Riccati 
equations and Sturm-Liouville problems has been firmly known at least since 
more than half a century [4] and more recently led to supersymmetric (SUSY) 
quantum mechanics [5J O [7] . 

In the present work we will discuss in some detail the periodic Sturm- 
Liouville (PSL) problems generated by particular solutions of the very simple 
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Riccati equations of constant coefficients [5] 

^ + R 2 + k 2 = (1) 
ax 

and 

^ - 2S<S> + $ 2 + k 2 + S 2 = , (2) 

which is closely related to the first one. We obtain the exact PSL solutions in the 
Bloch form. In general, when available, the analytic solutions of a PSL equation 
expressed as Bloch functions allow one to analyze the band structure of the 
problem. Of course, this is not always possible. In many cases, even if the exact 
solutions are known their Bloch form is hard to disentangle. As a consequence, 
other well-established approaches are frequently used in order to analyze the 
band structure of the spectrum: either by means of the Hill discriminant (or 
Lyapunov function) [3 or using the band structure parameter introduced by 
James [S]. In general, for this goal, two linearly independent solutions are 
required for all values of the spectral parameter. Hence, it is essential to have 
an as simple form as possible for the sought solutions in terms of the spectral 
parameter. Therefore we show that a convenient study of the band structure of 
the spectrum can be achieved using a representation for solutions in the form 
of power series with respect to the spectral parameter. 

The outline of the paper is the following. The next section is devoted to the 
three main points: (i) the SUSY construction of the periodic Sturm-Liouville 
problems of Schrodinger type starting from the aforementioned Riccati equation, 
(ii) solving the resulting Schrodinger equation and (iii) studying some of its prop- 
erties. Next, in Section 3, we provide a numerical approach of the same problems 
using the self-matching method of H.M. James [9] for the Kravchenko repre- 
sentation of the solutions in terms of spectral parameter power series (SPPS) 
[TOlITT]. A small Conclusion section ends up the paper. 

2 Riccati-associated PSL equations 
2.1 Case I: Riccati equation (OQ) 

In the Riccati equation fl} we introduce R(x) = 4- leading to the second order 
linear differential equation 

- f"(x) + vif{x) = , (3) 

where v\ — — fc§. 

The linearly independent solutions of ([3]) are obviously 

fi{x) = cosfcox , f2(x) — smkox . 

Using fi one gets for the Riccati equation the solution Rf 1 (x) — —kotankox. 
In what follows we employ fi to define Rf 1 (x), equally well one can use 
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leading to Rf 2 (x) — k cotk x, however this brings in only minimal changes in 
the whole of the mathematical apparatus that follows and hereby we will deal 
only with the first choice. 

Since v\ can be calculated by employing the equation 

the supersymmetric partner potential will be given by 

dR 

v 2 {x) = + R 2 h = fc 2 (l + 2tan 2 fc :c) 

and the supersymmetric partner equation of equation ([3]) is given by 

-g"(x)+v 2 (x)g(x)=0 . (4) 
The supersymmetric partner linear independent solutions are 



9i( x ) = ZZ~i — > 92{x) 



cos kox fcg cos kox 



k x 1 . 

— - — h — sin zkqx 



Considering now the spectral issue for these two periodic potentials, the 
constant one V\ and the singular one vi(x): 

-f"(x)-k 2 f(x)=K 2 f(x), 
-g"{x) + fcg(l + 2tan 2 fc a;).g(a;) = K 2 g(x), 

we can get the K 2 spectrum of the / problem from the Bloch solutions e ±W k o+ K2x 
that provide the quasimomentum P(K 2 ) — \/k\ + K 2 . On the other hand, 

for the g problem, being the Darboux partner of the / problem, i.e., g{x) = 
(<E — Rfi( x ))f( x )i we have the following Bloch solutions with the same quasi- 
momentum P{K 2 ) 



As known [T2] , the allowed energy bands exist only for P G R leading to K 2 > 
— Kq, therefore there is only one forbidden zone covering the interval (— oo, — fcg). 

Using the following two Pauli matrices 

-i \ . / 1 

1 o ) aid 0, 

a single matrix equation for the two-component spinor W = I J~ ) can be 

V h J 

written down: 

o-yW'(x) + ia x R fl W(x) = . (5) 

The two components of the spinor W have a sort of hidden coupling through 
the same particular Riccati solution but otherwise they look independent one 
of the other. 



3 



2.2 Case II: Riccati equation ([2]) 

A simple way to perform a direct coupling of the two components / and g is 
to add two constant parameters as follows: K as a spectral parameter of the 
Dirac-like equations and a potential coupling parameter S 



dx 



+ ia x (R fl + S)]W c {x) = iKW c (x) , W c = 9 £ J . (6) 



For related mathematical procedures the reader is directed to a textbook of 
Lanczos |13j and to a paper of Nogami and Toyama [14j and references therein 
for similar supersymmetric structure of the Dirac equation in particle physics. 
Equation ([6]) is equivalent to the following system of coupled equations 

g' c (x) + (koR fl + S)g c (x) = Kf c (x) (7) 

f' c (x) ~ (k R fl + S)f c (x) = -Kg c (x). (8) 

This leads to: 

s + *) (sir ••)'•-*"*>■ < 9 » 

(10) 

where $i = Rf y + S = — fcotanfcoir + S is a particular solution of the Riccati 
equation In unfactorized form, ^ and (TTUj) turn into the following equations 

-/" (z) + (-*§ + 5 2 - 2SManfcox)/ c 0r) = Jf 2 /cW (11) 
-5c ( x ) + ( fc o + ~ 2Sk ta,nk x + 2fc 2 tan 2 k x)g c (x) = K 2 g c {x) . (12) 

The latter two equations define two new classes of parametric singular po- 
tentials. To see the changes with respect to the initial equations ([3]) and ([4]) 
for / and g, respectively, we write the previous system in terms of modified v 
functions v\ c (x) and v-i c {x): 

-K(x) + v 1>c (x)f c (x) = K 2 f c (x), -g'^x) + u 2>c (x)g c (x) = K 2 g c (x) , (13) 

where 

vi, c (x) = vi+Av(x), V2,c{x) = V2{x)+Av{x), Av(x) — S 2 —2Skot&nkoX . 

Note that all the solutions f c of (fTTj) for any value of K 2 are square integrable 
on any finite interval. This is due to Weyl's alternative theorem, see for example 
the book of Hellwig [T5], all singular points (n odd) of v\^ c (x) are limit 
circle points. Indeed, it is sufficient to prove square integrability of two linearly 
independent solutions for any fixed value of K 2 . Taking K 2 = we see that 
both linearly independent solutions 

f e~ 2Sx 

hi(x) — e Sx coskox and h2(x) = e ^cos/co^ / t~, — dx 

J cos 2 k x 
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are square integrable on any finite interval. 

The limit circle points of equation become limit point singularities 
of (fT2)) . This can be demonstrated considering the solution of P^|) for 

K 2 = 0. Thus, Weyl's alternative guarantees that exactly one square integrable 
solution of (|12p exists for K 2 with a nonvanishing imaginary part whereas for 
real K 2 such a solution cannot even exist. 

In order to solve the equations (fTTjl . (fT2|) we begin with the first equation in 
(fl~3|) . With the aid of its solution f c the solution of the second equation in (|13| 
can be simply obtained by applying the Darboux transformation 

g c (x)=f c (x)-$ 1 (x)f c (x). (14) 

Thus, we focus next on the f c equation which will be reduced to a hypergeo- 
metric equation. The reduction is done with the help of a procedure similar to 
the one described in [16] , 



2.3 Hypergeometric solutions 

Consider the simultaneous change of the independent variable x = \ (1 — * tan kox) 

2 COS l x and of the dependent variable y(x) = (x 2 ~x) 1 ^ 2 /c- Then the first equa- 
tion in (fT5|) takes the following Schrodinger-like form 



(15) 



where 



S 2 - K 2 + 2iSk - MSkpX 
4fc V(x-l) 2 



Finally, in order to bring (fT5"]) to the hypergeometric form, the following substi- 
tution can be used [PT] 



where 



y(x) = x p (x - i) q u( x ) , 



Pl,2 



1± 



! (k -iS) 2 +K 2 



91,2 



1 ± 



l(k + iS) 2 +K 2 
k 2 , 



Thus, one gets 

^d 2 U \ , d U 

xix - l )-rr + 2 Kpi,2 + qi,2)x - pi,2\ -y- 



dx 



dx 



2pi,29l,2 



K 2 



2k 



U = (16) 



Choosing the pair pi and q\ , we obtain the following solutions 

Ui = a Fi(pi +gi-l,pi+<7i; 2 Pl ; X ) , U 2 = x^ 2pi 2 F 1 {- Pl +q 1 +l, -pi+qr, 2-2 Pl ; x) ■ 

(17) 



5 



Using properties of hypergeometric functions (the change of variable \ ~ * v^i)' 
the f c linearly independent solutions can be written in the form 

f cA (x) = (-l)<K-k e -i{^-^o* 2Fl ( pi + qi -l :Pl - qi; 2 Pl ; - e ~ 2 ^) , 

(18) 

fed*) = (-lyi-ie^- 1 ^ * 2 F 1 (- Pl + qi ,- Pl - qi + l;2- 2 Pl ;-e~ 2ik ° x ) . 

(19) 

Application of the Darboux transformation ([M]) to the solutions f c ,i(%) and 
f c ,2{x) leads to the following solutions g c ,i{%) and g c ,2(x) of the second equation 
in' US]) 

9c,l(x) = (-l)«i-5 e -<2 P1 -l)fco* X 

x [(fc tanfc x- > S'-ifc (2p 1 -l)) 2 F 1 (p 1 + gi -l,p 1 - gi ;2p 1 ;-e- 2jfcoa; ) + 

+ aF (pi + _ + 1; 2pi + 1; _ e - 2ifcoa:) 

pi 

(20) 

<7c, 2 (a;) - (-lJfi-ie'Cto- 1 )*"^ 



[(fc tanfco^ - 5 + ifc (2pi - 1) ) 2*1 + Ql , - Vx - q x + 1; 2 - 2 Pl ; _ e - 2 *>*) + 

( Pl + g x -l)( P i-gi) . fcoe _ Mfeoa ( + + _ + 3 _ _ e _ 2ife0S) 

1 - Pi 

(21) 



Though the Darboux transformation is applied to regular solutions / C) i and / C) 2 
the singularity of the superpotential $ implies the singularity of g cA and g Ci 2- 

The solutions (JTHJ), (HU) and (j2T))) . (f2"Tj) are quasiperiodic or Bloch functions 
with the quasimomentum 



P S (K 2 ) = (2pi - l)fe = V(fco - iS) 2 + K 2 , 

which defines the Brillouin zone as follows: Re(Ps) £ [— fco,^o] ■ The allowed 
energies exist only for Pg G R [12]. This condition holds when (ko — iS) 2 +K 2 £ 
R + . To specify the spectrum we should make some additional considerations. 
Limiting ourselves to the real values of the spectral parameter K 2 we have that 
for ko £ TL, S should be of the form 5 = is, where s £ R. Moreover, K 2 must 
satisfy the inequality K 2 > —(ko + s) 2 or equivalently K 2 £ (—(ko + s) 2 ,oo]. 
It is worth mentioning that for s = we get the spectrum of the uncoupled 
potentials v\ and v 2 . Another apparently possible case: S £ R, ko = vy, where 
7 £ R that also leads to a real Ps is not meaningful for the analysis of issues 
related to the quasimomentum since the potentials Ui tC and v 2%c are not periodic 
any more. Thus the spectrum of Ijlip and (|12p is real only for purely imaginary 
values of the parameter S. 
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Notice that periodic potentials with purely imaginary coupling constants 
have been considered for the case of Mathieu equation [18] and have applica- 
tions to the alternating flow of electromagnetic fields along conducting elliptic 
cylinders. 

{2nkn 
, ,n = 0,±1,±2, • • • . The 

(Zn + 1 )kq 

solutions / c> i and f Cj2 are periodic when 

P s = 2nk Q =>■ K 2 = (2nk Q ) 2 - (k - iS) 2 (22) 

and antiperiodic for 

P s = (2n + l)k =^K 2 = ((2n + l)fc ) 2 - (k - iS) 2 . (23) 

For these special values of K 2 corresponding to the band edges two pairs of 
Bloch solutions f c ,i(x, K 2 ) 7 f c _ 2 (x,K 2 ) and g c ,i(x, K 2 ), g c ,2(x,K 2 ) degenerate 
to single solutions f c .i(x,K 2 ) = f Ci2 (x,K 2 ) and g c ^(x,K 2 ) = g c , 2 (x, K 2 ). In- 
deed, the solutions (fTS)) and (fT^) in the periodic case Ps = 2nk^ have the form 



= (-iy^-h- 2ink °* 2 Fx (n-\ + qi(n),n + 1 - qi (n); 1 + 2n; _ e ~ 2 *^ , 

fc%(x) = (-l)^ m )--2e 2lmkax 2 Fi (-m - ~ + qi (m), -m + ~ - gi (m); 1 - 2m; _ e - 2 *^ , 
where qi{n) — \ + \Jn 2 + i-^. 

It is clear that in order to have the same value of Ps in both functions we 
should take m = — n, thus the solutions are periodic and /^(i) = f£ %{x). 



In the antiperiodic case Ps = (2n+l)ko we have q\ (n) = i + J (n + ^) 2 + i-j^ 
and 

/cT» = (-l)^-h-^ 2n+1 ^ x a Fi (n + <ft(n),n + 1 - ?i(n);2n + 2; - e - 2 ^) , 

/ c T r ( x ) = (-l)« l(m) -Se i < am+1) *° a 2 i ; i (-m-l + qi {m),-m- qi {m)--2m--er 2lkaX ) 

Taking m = — (»+ 1) we obtain that /"i 6r = /" ^ are antiperiodic. Note 
that in both cases qi(n) = qi(m). Analogously we obtain 

g p c er {x) = (-iyiW-h e - 2inkoX \ (k teuik x-S - 2mk ) x 

x 2 Fx (n - i + qi(n),n + ^ - q 1 (n); 1 + 2n; _ e - 2 *^ + 

S_ e -2 lkoX ^ + 1 + gi(n)j n + 3 _ ?i(n) . 2 + 2n . _ e -2i fc0 x^ 



n , 2 
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and 

g a c per {x) = (_l)9i(n)-| e -i(2n+l)feox \ ^ tan _ £< _ ( 2n + l)ifc ) X 

x 2 F 1 (n + q 1 {n),n + l- ?1 (n); 2n + 2; - e - 2 *>*) + 

S e^' * 2 Fx (n + 1 + gi(n), n + 2 - ?1 (n); 2n + 3; -e- 2 **"*) 



n + 1 



The value Ps = and consequently i^g = — (fco + s) 2 give us the following 
periodic nodeless solution of (fTTj) which we denote by fo(x) 



f (x) = (-iyV^ 2 F! ^ 

This eigcnfunction will be used later on 



3 An efficient numerical approach for the energy 
band structure 

The potentials i^i iC and V2, c are periodic functions of period A = and have 
singularities at the points j^-. Following James [9] we choose the first period 
as [0, A] and call it the zeroth cell, the second period [A, 2A] as the first cell, 
and so forth. Following the fundamental procedure of James [3], we construct 
the so-called self-matching solutions of the SUSY-related equations ((TTj) and 
(fT2"|) for the zeroth cell which allows us to build the Bloch solutions on the 
entire range of x. We proceed further by writing the SPPS representation of 
the associated Hill discriminants which allows us to describe the spectrum of 
the SUSY-related equations (fTTj) and (fT2")) in a simple way. For doing this we 
choose to use a numerically calculated solution instead of the exact one given 
in the preceding section. First, because the SPPS approach is clearly more 
universal and can be applied in situations when the exact solution is unavailable. 
Second, and more important for this work is that all the following constructions 
imply the computation of solutions for a large set of different values of the 
spectral parameter K 2 , while the exact solutions involving the hypergeometric 
functions have been proved considerably less practical than the approximate 
solutions obtained below. The SPPS method allows one to construct a solution 
in the form of a power series with respect to the parameter K 2 which is ideally 
suited for our purposes. Compared to the use of exact solutions it gives us the 
possibility to calculate a solution for different values of K 2 in a more efficient 
way. From Figs. (1) and (2) one can assess the excellent agreement displayed 
by the zeroth-cell solutions obtained using both methods for a given set of the 
spectral parameter K 2 . 

3.1 Self-matching cell solutions 

Now we begin with the equation for f c on the zeroth cell x G [0, A] in order to 
construct the so-called self- matching pair of independent cell solutions F±[p_]. 
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To obtain these special solutions it is necessary first to have two linearly inde- 
pendent solutions /i and fi satisfying the following initial conditions 

/i(0,iT 2 ) = l, / 2 (0,^ 2 ) = 0, (25) 
f[(0 1 K 2 ) = 0, f' 2 {Q,K 2 ) = l. 

The method of spectral parameter power series (SPPS) [TQl [11] gives these 
solutions in explicit form as follows. Let /o be a particular solution of — /q + 
i^i/o = Kq/o such that f S C 2 (0,A) together with are bounded on [0, A]. 
The general solution f c has the form f c = G\fi + C2 J 2 , with 

(26) 

/a(i)=-/o(0)/o(a!)S l (i) ) 

where E andEi are the spectral parameter power series E (x) — J2°^=o X^ 2n \x){K 2 
K$) n , Ei (a;) = E^i ^"^^X-^-^o)™ -1 with the coefficients X, X given 
by the following recursive relations 



'^^^(O/oCO* foranoddn 
A (n) (x) = { (28) 



j^(n-i)^_^_ for an even n 



fj* X (n-i)^_d£_ foranodd) 
X^ n \x) = \ ' (29) 

[f* X( n - l \Z)fg(Z)dt foranevenn . 

The solution fo is given by and corresponds to the particular value of 
AT, 2 = — (ko + s) which represents a band edge. 

One can check by a straightforward calculation that the solutions f± and 
/2 fulfill the initial conditions ([23]) . In Fig. 1 we plot the series solutions /1 
and f% evaluated by (f2"6"]l by a solid line and the markers represent the exact 
solutions calculated taking the appropriate linear combinations of / Cj i and / C) 2 
given respectively by {15]) and ([TO]) in order to fulfill the same initial conditions. 
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Fig. 1. The solid lines represent the series solutions /i and ji evaluated by (|26[1 and 
the markers represent the exact solutions (|18p and (|19|l . The values of the parameters 
S and ko are O.li and 1, respectively. 

Denoting now 

h(A,K 2 ) = ail (K 2 ), f 2 (A,K 2 ) = a 12 (K 2 ), 
/{(A, K 2 ) = a 21 (K 2 ), K(A,K 2 ) = a 22 (K 2 ), 

one can write the self-matching solutions F± in the form 

F±(x, K 2 ) = h(x, K 2 ) + a±f 2 (x, K 2 ), 

where a± solve the equation ai 2 a± + (an — a 22 )a± — a 2 \ = [9]. 

To obtain the self-matching solutions to the equation (fT2|) we first con- 
struct the solutions g\ and g 2 which satisfy the initial conditions <7i(0, K 2 ) = 
g' 2 (0,K 2 ) = 1 and g 2 {0,K 2 ) = g[(0,K 2 ) = 0. For this taking the following 
linear combinations 

hi*) = ^ ( s h( x ) + ( s2 + r2 ) M x )) and U x ) = (fx(x) + Sf 2 (x)) 

and applying the Darboux transformation (fT4]) to them gives 

gi (x) = f[(x) - and g 2 (x) = f' 2 {x) - ®x{x)f 2 {x). (30) 

Illustrative plots of the latter solutions are displayed in Fig. (2) in solid lines, 
while the markers correspond to the exact formulas ([2TJ)l and (f2"TT) . Notice also 
that all the singularities of the above solutions are contained in the Darboux 
transformation function <5>\{x) = S — feptanfeo^- 
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Fig. 2. Series and exact solutions gi and gi evaluated by (|30[l and (|20l) and (|21[) . re- 
spectively, for the same parameters as in Fig. 1. 

By analogy with the /-case we denote 



where (3± are roots of the equation 6i2/3± + (&n — b 2 2)f3± — 621 = 0. 
3.2 Bloch solutions 

We are now in a position to write down the Bloch (quasi-periodic) solutions to 
the equation (fTTj) through the whole range of x divided as follows nA < x -< 
(ra+l)A for n = 0,±1,±2,- • • 



f±(x,K 2 )=rlF ± (x-nA 1 K 2 ) = rl [fi{x - nA,K 2 ) + a±f 2 (x - nA,K 2 )] . 

(31) 



The Bloch factors r± are a measure of the rate of increase (or decrease) in 
magnitude of the self-matching solutions F±(x, K 2 ) when one goes from the left 
end of the cell to the right one, i.e., 



gx(A,K 2 )=b n (K 2 ), g 2 (A,K 2 ) 
g' 1 (A,K 2 )=b 21 (K 2 ) 1 g' 2 (A,K 2 ) 



bMK 2 ), 
b 22 (K 2 ), 



and the self-matching solutions G± have the form 



G±(x, K 2 ) = 9l (x, K 2 ) + (3±g 2 (x, K 2 ), 



r±{K 2 ) 



F ± (A,K 2 ) 
F±(0,K 2 ) 



The values of r± can be also written as 



r ± (K 2 ) = - (D f {K 2 ) T yfD*(K*)-4) 
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where Df(K 2 ) denotes Hill's discriminant (also known as Lyapunov function) 
associated with (fTTj) [3] and Df(K 2 ) = an + a 22 - For equation (fl~2j) the Hill 
discriminant is given by D g (K 2 ) = bu + b 22 . Using the relations ([50)1 between 
the solutions g\, g 2 and /1, /2 the identity Dj{K 2 ) = D g (K 2 ) can be easily 
obtained. This means that the Bloch factors for the quasi-periodic solutions 
to equation (fT2")l are the same as for (fTTj) . Thus, for a numerable set of cells, 
nA < x -< (n + 1)A for n = 0, 1, 2, • • • , one can write 

g±(x,K 2 ) =rlG ± {x-nk,K 2 ) = r n ± [ 9l (x - nA, K 2 ) + f3 ± g 2 (x - nA, K 2 )] . 

(32) 

In Fig. 3 the Bloch solutions (|3"Tj) and (f3"2"| are plotted with solid and dotted 
lines respectively. 




Fig. 3. The blue lines are the solutions f± given by (|31[1 and the red lines represent 
the solutions g± given by (|32[1 for the same values of the parameters as previously and 
for K 2 = 0.25 (solid lines), 2.25 (dotted lines), and 6.25 (dashed lines). 



3.3 A power series representation for the Hill discriminant 

The Hill discriminant (Lyapunov function) allows one to describe the spectrum 
of periodic differential equations. Namely, the spectrum of (fTT|) and (fT2|) is given 
by the following set [2] {K 2 : Df{K 2 ) g R and |Z?/(.?r 2 )| < 2}. The expression 
for Df(K 2 ) = &n+a22 can be written in a simple explicit form. For this we write 
an and 022 in form of a spectral parameter power series using (|26p and taking 
into account that £^i(x) = - r j^- ) Z (x), where £ (z) = J2n=o x(2n) (x){K 2 - 

K 2 ) n : 

a n = M^S (A)+/^0)/o(A)S 1 (A) and a 22 = -/o(0)/^A)S 1 (A)+^^E 

Since /o(a;) is a A-periodic function: /o(0) = /o(A). Finally, writing the ex- 
plicit expressions for Eq(A) and So (A) we obtain a representation for the Hill 
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discriminant associated with (fTT|) and (fT2)) 

oo 

Df(K 2 ) = Z?,,^ 2 ) = (^ (2n) (A) + X( 2 ")(A)) {K 2 - K 2 ) n . (33) 

Thus, only one particular nodeless periodic solution fo(x) of (fTTj) is needed for 
the construction of the Hill discriminant Df(K 2 ). There are other known series 
representations of the Hill discriminant, see [31 [TS] . Nevertheless none of them 
allows one to represent it as a spectral parameter power series which is extremely 
useful for calculations involving different values of the spectral parameter. 

Figure 4 shows the plot of Df(K 2 ) which we evaluate in two ways. With the 
solid line we plot the function Df(K 2 ) obtained by means of the SPPS solutions 
given by (|33| and the markers correspond to Df(K 2 ) obtained with the exact 
solutions (HHJ) and JTSJ). From (O the advantage of the SPPS method for 
calculating the Hill discriminant and hence the corresponding spectrum can be 
assessed. Using the SPPS the calculation of the value of the Hill discriminant 
for every value of its argument reduces to a simple substitution of the value 
of K 2 into an easily evaluated expression (j33|) . The values of K 2 for which 
Df(K 2 ) — ±2 correspond to the band edges of the spectrum. Notice that 
D f (K 2 ) = ±2 is in accordance with (J22J) and l[2"3"j). namely D f (K 2 ) = 2 exactly 
for K 2 = (k -iS) 2 -(2nk ) 2 and D f (K 2 ) = -2 when K 2 = (k - iS) 2 - ((2n + 
l)fc ) 2 . 

When D(K 2 ) ^ ±2, the general solutions of (fTTj) and (fT2]) have the form 
(understanding that it refers to the whole x axis henceforth) 

f c (x) = C+f+(x) + C-f-(x) and g c (x) = C+g+{x) + C^g.{x) 

For the values of K 2 giving Df {K 2 ) = ±2, two pairs of independent solutions 
f + (x,K 2 ), f^.(x,K 2 ) and g + (x,K 2 ), g-(x,K 2 ) reduce to a pair of a single 
solutions f + (x,K 2 ) = f-(x,K 2 ) and g + (x,K 2 ) = g-(x,K 2 ) , but there is a 
definite prescription for the construction of an independent second solution [9] . 
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Fig. 4. The Hill discriminant Df(K 2 ) evaluated using SPPS solutions ((26) (solid line) 
and by the exact solutions (|18[1 and (|19|) ( markers) for S = {0.1i,0.7i, 1.5i}. 

4 Conclusions 

We have used one particular solution of simple Riccati equations of constant 
parameters to build the corresponding supersymmetric partner Sturm-Liouville 
equations. The latter equations are solved analytically in terms of hypergeo- 
metric functions. Furthermore we worked with Kravchenko's spectral parame- 
ter power series solutions that are better suited from the algorithmic (numer- 
ical) standpoint and allows an easy implementation of the old self-matching 
procedure of H.M. James [9] for solving periodic Sturm-Liouville problems of 
Schrodinger type in terms of Bloch solutions. We also obtain an effective power 
series representation of the Hill discriminant in terms of the Kravchenko series. 
The mathematical procedure expounded in this paper can be applied to more 
general periodic SL equations that abound in the area of nanostructured ma- 
terials and in the form of periodic Hclmholtz equations in photonics. Other 
applications can be foreseen in the areas of chirp technology, see for example 
Refs. HH and [5T]. 

The first author would like to thank CONACyT for a postdoctoral fellowship 
allowing her to work in IPICyT. 
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